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J BESSEL FUNCTION (FORTRAN) 
DECUS Program Library Write-up DEC US No. 8-96 



C; This i^KOGKA.vl CO:^iPUr£b Th F: J RFSSEL FUNJCTION 

C; FOR A GIVFN AkGU'^jENT AND OKDFH 

C; FfVKOH CODE key: 

C; IER = C1 -NO FRROR 

C; IF.R=1.. ORDER I !S NEGATIVE 

C; IER = S ARGUMENT IS >\EGAT1VE OR ZERO 

Ci IER = 3 REO'D ACCURACY NOT OorAI.\ED 

Ci IFR = 4 RANGE OF N* wRT. X NOT CORRECT 

Si^2; CONTINUE 

TYRE 100 
100; FORviAT (/***EMTER ORDER OF J BFbSFL FUMCTION DESIRED"*/) 

ACCEPT 200 *i\ 
200; FORMA TCI) 

TYRE 300 
300; FORviATC'ENTER ARGUMENT OF J dES.SEL FUNCTION DESIRED"*/) 

ACCEPT 400 *X 
4005 FORNIATCE) 

TYPE 500 
S00; FOR'MATC "ENTER ACCURACY DESIRED"*/) 

ACCEPT 400*1) 

BJ=0. 

IFCN) 10*20*20 
105 IER=1 

GO TO 77 
20 5 IF(X)30*30*31 

30 5 IER=2 

GO TO 77 
315 IF(X-1 5.)32*32*34 

32 5 NTEST=20.+10.*X-X**2. /3. 

GO TO 36 
3 45 NTEST=90.+X/2. 

36; IF(N-NTEST)^(^*3H*38 

3H5 IER=4 

GO TO 77 
40; IER=0 

Nl =N+1 

8PREV=.0 

IFCX-b. ) 53* 60*60 
50 5 v,A = X + 6. 

GO TO 70 
60; vjAsl . 4*X + 60 . /X 

lFiX=X 
705 Mrt = \i+IFlX/4 + 2 



\]2; 

111; 
1 13; 



110; 

1 30 ; 
1 30 ; 



1 40 5 

1 b0; 



1 m; 



1 70; 

1 80 ; 



191; 
1 9P.i 

193; 
1 9 4; 
1 95J 

1 90 5 

77; 
10S; 



I F ( virt - v)A ) n 1 , 11 >> , 11 2 

GO TO 1 1 3 

C ON T IN) UK 

•v!.v!AX = N'rFSr 

1)0 190 !^=MZF:f^0*'^i^^AX>3 

F^il =1 .0R-PH 

K^=ALPHA=0. 

IKC'^-C v;/0)*0) 1P0^ 1 10;. 120 

JT = -1 

GO TO 1 30 

JT=1 

:V) 2 = M - 2 

00 160 K=l,"^?> 
!v)K = vi-K 
FLOAT = ^^K 

r?viK=2. + FL0Ar*F'^<l /X-F^1 

F'VjrrFVil 

1 F ( ^K - \'- 1 ) 1 50 > 1 40 , 1 50 

jT=-J1 
b=l+JT 

ALPHA = ALPHA + rS.-^K + S 

CO^TINUF 

R^^K = 2.*Fivil /X-F^ 

IF(N) 180> 1 7(1* 180 

ALPHA = ALPHA+Bv,K 

BJ=BJ/ALPHA 

t.ltJ=^3J-HPf^FV 

IFC UU) 1 91 > 191 > 192 

uu=-uu 

IFC VV) 1 93* 193* 194 

v\;=-vv 

IF(UU-V\/)77*77* 195 

rfPREV=Bj 

CO.MriN'UF 

1 Fk = 3 

TYPE 102*BJ*1P:H 

FORV)AT(F*" 1FP = "*1*/) 

GO TO 502 

END 



